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Abstract 



The influence of molecular vibrations on dissociative adsorption is stud- 
ied by six-dimensional quantum dynamical calculations. For the system H2 at 
Pd(100), which possesses non-activated pathways, it is shown that large vi- 
brational effects exist and that they are not due to a strongly curved reaction 
path and a late dissociation-hindering minimum barrier, as was previously 
assumed. Instead, they are caused by the lowering of the H-H vibrational fre- 
quency during the dissociation and the multi-dimensionality of the potential 
energy surface. Still there are quantitative discrepancies between theory and 
experiment identified. 
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I. INTRODUCTION 



In recent years the dissociative adsorption and associative desorption of hydrogen on 
metal surfaces has been the subject of many experimental and theoretical investigations 
(see, e.g., Refs. [|l],0|| and references therein). The studies were motivated by the fact that 
for these relatively simple systems the process of bond breaking and bond making during 
dissociative adsorption can be analyzed in greater detail. 

Breaking a molecular bond is obviously coupled to the interatomic distance and molec- 
ular vibrations. Therefore a large number of studies have addressed the influence of those 
vibrations on the dissociative adsorption and associative desorption, especially in the bench- 
mark system H 2 /Cu @,§|,P,0,PJ1] • In desorption studies strong vibrational heating of 
the hydrogen molecules was found, i.e., vibrational population ratios were much greater than 
expected for thermal equilibrium with the temperature of the substrate at which desorption 
occurs. According to the principle of microscopic reversibility this implies that the proba- 
bility of dissociative adsorption should be enhanced for vibrational excited molecules, which 
indeed has been confirmed experimentally ||0f^||. These vibrational effects were usually 



discussed [|IT)| , |lT||i2| , |T3| within the context of potential energy surfaces (PESs) which depend 
only on two coordinates, namely the center of mass distance of the molecule from the sur- 
face, Z, and the H-H interatomic distance, d. Obviously, the motion of the two atoms of a 
diatomic molecule is governed by six coordinates, not just Z and d, but the dependence on 
the other (neglected) four coordinates was felt to be of minor importance. It was generally 
accepted that the PES in the considered two-dimensional space (a so-called "elbow plot") 
should exhibit a strongly curved reaction path and a so-called late barrier towards disso- 
ciative adsorption, i.e. a barrier after the curved region of the PES, close to the surface, 
in order to account for strong vibrational effects in dissociative adsorption and associative 
desorption. 

Interestingly, strong vibrational heating was also found for hydrogen molecules desorbing 
from Pd (100) although for this substrate, in contrast to Cu, the dissociative adsorption 
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is non-activated po|,p!o[pT| . Nevertheless, within the spirit of the above discussion for the 
H 2 /Cu system, Brenig et al. [|18|,[nj reproduced the vibrational heating for H 2 desorption 
from Pd (100) by quantum dynamical calculations using two-dimensional model potentials 
with a minimum barrier of 200 meV. Since this is in conflict with the fact that adsorption of 
H 2 in non-activated, Darling and Holloway |2T)| , |2"T| questioned the validity of this theoretical 



work and argued that one has to take into account a distribution of barrier heights which 
requires higher-dimensional calculations. In their model calculations, which were still two- 
dimensional, they basically showed that it is not possible to reproduce vibrational heating 
in desorption with a barrier-less two-dimensional elbow potential [pCfl . 

Recently we enhanced the computational approach for scattering of molecules at sur- 
faces |L9"l , p2"! and are now able to investigate reactions of diatomic molecules on surfaces 
with all six degrees of freedom of the molecule being treated quantum dynamically |f23| . 



Such a study was performed for the sticking probability of H 2 at Pd (100) employing a 
high-dimensional PES derived from first-principles calculations ||17|| . In the present paper 
we use the same PES and the same six-dimensional (6-D) quantum-dynamical method. The 
PES for the interaction of H 2 and Pd (100) possesses non-activated pathway, but also (in 
fact in majority) activated pathways []17|| . It was found that H 2 molecules impinging with 
small energies are efficiently steered along the non-activated pathways towards dissociatice 



adsorption [23]. However, with increasing energy the steering effect gets less effective and the 



majority of molecules proceeds via pathways with energy barriers [23]. The steering mech- 



anism is operative for the translational degree of freedom as well as for the rotations p4| 



In the meantime two predictions of the quantum dynamical calculations for H 2 /Pd(100), 
namely the strong rotational hindering of the dissociation at low kinetic energy and 



the orientation dependence of adsorption and desorption |23" |, have been confirmed experi- 
mentally |2|J2] 



In order to investigate the effects of molecular vibrations on the sticking probability we 
had to extend the calculations P3[ to higher energies. To describe the dynamics of initially 



vibrating molecules properly we have increased the number of vibrational eigenfunctions in 
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FIG. 1. Contour plot of the PES along a two-dimensional cut through the six-dimensional 
coordinate space of H2/Pd(100). The inset shows the orientation of the molecular axis and the 
lateral H2 center-of-mass coordinates. The coordinates in the figure are the H2 center-of-mass 
distance from the surface Z and the H-H interatomic distance d. The dashed line is the optimum 
reaction path. Energies are in eV per H2 molecule. The contour spacing is 0.1 eV. 

the expansion of the wave function of the hydrogen nuclei: In the present study up to 25,200 
channels per total energy are taken into account compared to 21,000 channels which were 
considered previously P3| . We will show that in spite of the absence of a minimum barrier 
towards dissociative adsorption and a strongly curved reaction path there are still substantial 
vibrational effects in the adsorption and desorption of H 2 /Pd(100). We will demonstrate 
that they are caused by a strong lowering of the H-H vibrational frequency during the 
adsorption and by the multi-dimensionality of the PES relevant for the dissociation process. 



However, some discrepancies to the experimental results of refs. [14.18] remain. 



II. RESULTS AND DISCUSSION 



Figure [T] shows a cut through our PES of H 2 /Pd(100), where the most favorable path 
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towards dissociative adsorption is marked by the dashed line. As discussed above, for this 
path there is no energy barrier hindering dissociation, i.e., the adsorption is non-activated. 
The curvature of the optimum reaction pathway (see dashed line in Fig. [I]) is relatively mod- 
erate compared to the curvature of reaction pathways of previously assumed or guessed PESs 
which had been used in earlier low-dimensional studies pTD|JTT|JT^JT^JT5|JT9|j^D[ ] . The detailed 
total-energy calculations showed that the PES is strongly anisotropic and corrugated so 
that besides non-activated paths the majority of pathways towards dissociative adsorption 
has in fact energy barriers with a rather broad distribution of heights and positions. The 
barrier height depends on the molecular orientation and impact site in the surface unit cell. 

Figure ^| presents results for the sticking probability as a function of the kinetic energy of 
a H 2 beam under normal incidence. Quantum mechanically determined sticking probabilities 
for hydrogen at surfaces with an attractive well exhibit an oscillatory structure as a function 



of the incident energy p3|j2^j28|j29|| , reflecting the opening of new scattering channels and 
resonances as a l so observed in, e.g, LEED [jSOfl . Such oscillations have not been 

observed yet. In the experiments the molecular beams are not strictly mono-energetic but 
have a certain energetic broadening. For the calculations of Fig. ^| we assumed an energy 
width AEi/Ei = 0.2, typical for experiments [T5|]. As a consequence, the quantum dynamical 
oscillations are smoothed out. 

The solid curves correspond to 6-D calculations for H 2 molecules initially in the vibra- 
tional ground and first excited state, respectively. The z/, = curve shows the characteristic 
initial decrease of the sticking probability with increasing energy (for < 0.1 eV) which is 
due to the decreasing importance of the steering effect with increasing energy [p3fl . Also for 
molecules initially in the first excited vibrational state a corresponding behavior is found, 
but for these the steering effect is strong only for very small energies (Ei < 50 meV). For 
kinetic energies higher than « 50 meV the Vi = 1 molecules experience an increasing sticking 
probability which is significantly larger than for z/j = molecules. 

The effect of the initial vibrational state can be quantified by the vibrational efficacy 
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FIG. 2. Sticking probability versus kinetic energy for a H2 beam under normal incidence on a 
Pd (100) surface. The molecules are initially in the rotational ground state ji = 0, and in their 
initial vibrational states are Vi = (lower curves) and V{ = 1 (upper curves). The solid lines show 
the results of the 6-D calculations and the dashed lines are five-dimensional calculations where the 
vibrational degree of freedom is approximated by an adiabatic treatment (see text). 



where e(S) is the kinetic energy required to obtain the sticking probability S, so that S vi b(5') 
is the separation of the sticking curves for a certain sticking probability divided by the gas- 
phase vibrational quantum tkv v n, = 516 meV. In Fig. |]we have marked S = 0.7, which leads 
to a value of the vibrational efficacy of S vi b(5' = 0.7) = 0.75. This means that 75% of the 
vibrational energy is effective in promoting the dissociative adsorption, a value even higher 
than in the H 2 /Cu(lll) system ||. 

To clarify why and how molecular vibrations effect the sticking probability we performed 
calculations for a reduced coordinate space, namely allowing only for five degrees of freedom 
for the two hydrogen atoms. This reduction was achieved by keeping the molecules in their 
initial vibrational state. Although the vibrational state is kept fixed, the energy of the 
vibration, which is determined by the strength of the H-H interaction, will change along the 
scattering pathway. As can be seen in Fig. these five-dimensional results are very close to 
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FIG. 3. H-H vibrational frequency Hlv(s), potential Vq(s) and "vibrationally adiabatic poten- 
tials" (see Eq. 2) along the reaction path of Fig. Q. s = corresponds to the point of maximum 
curvature of the reaction path. 

the 6-D results. This reflects two facts. First, the molecular vibrational state is a sufficiently 
good quantum number and is almost conserved during the scattering, i.e., the probability 
for transitions between different vibrational states during the scattering event is rather low. 
And second, the curvature of the reaction path of the H 2 /Pd(100) PES is not crucial for the 
vibrational effects in this system because in the 5-D calculations no curvature is present in 
the Hamiltonian. 

As the next step we will analyze the "vibrationally adiabatic potentials" which are defined 

by 

KkOO = v (s) + (huj(s) - ^ vib ) („ 4 + i) , (2) 

where s is the coordinate along the reaction path (see Fig. |l|). The vibrationally adiabatic 
potential is the relevant potential for the H2 molecule moving on the PES in a particular 
vibrational state taking the change of the vibrational frequency along the path into account. 
In Fig. |3| we have plotted the vibrational frequency together with the potential and the 
vibrational adiabatic potentials for z/j = and = 1 along the reaction path coordinate s 
of Fig. [l]. At s = 0, the point of maximum curvature in Fig. [l|, the vibrational frequency is 
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strongly reduced from its gas phase value of hu = 516 meV to about 150 meV. This leads to 
a lowering of the vibrationally adiabatic potential by 183 meV for — and by 549 meV for 
Vi = \. Such a lowering does not only occur for the most favorable adsorption path, but also 
for other non-activated and activated pathways, i.e., for other impact sites in the surface 
unit cell and for other molecular orientations. This can be demonstrated by the integrated 
barrier distribution for the potential Vq and the vibrationally adiabatic potentials 



In Eq. H|, 6 and <p are the polar and azimuthal orientation of the molecule, and X and Y 
are the lateral coordinates of the hydrogen center-of-mass. A is the area of the surface unit 
cell. Each quadruple defines a cut through the six- dimensional space (see Fig. [I] for one 
example), and is the minimum energy barrier along such a cut. The function O is the 
Heavyside step function. The quantity Pb(E), which is plotted in Fig. f|, is the fraction of 
the configuration space, for which the barrier towards dissociative adsorption is less than E; 
it corresponds to the sticking probability in the classical sudden approximation, which forms 
the basic approximation behind the so-called "hole model" |3"I|| . Actually, the comparison 
of Figs. and f| reveals that the hole model gives a satisfactory description of the sticking 
probabilities at high kinetic energies above 0.3 eV, whereas it is strongly at variance with 
full dynamical calculations at low kinetic energies, where the steering effect is operative. 

Figure |] demonstrates that the barrier distribution is lowered due to the decrease of 
the vibrational frequency by about 180 meV for molecules in the vibrational ground state 
compared to the potential Vq(s) and by further ^400 meV for molecules in the first excited 
vibrational state. Because of the lowered potential the vibrationally excited molecules are 
accelerated more strongly towards the surface, they become so fast that the steering mecha- 
nism is suppressed. For that reason the initial decrease of the sticking probability for z/j = 1 
is limited to low energies (Ei < 50 meV, see Fig. 0). The difference in the barrier heights of 
about 400 meV is reflected by the energetic shift between the sticking curves for z^ = and 
z/j = 1 in Fig. ^|for sticking probabilities larger than ^0.7. Hence it is the strong decrease of 
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FIG. 4. Integrated barrier distribution using the potential energy Vq and "vibrationally adia- 
batic potentials" . 

the H-H vibrational frequency during the dissociation which causes the vibrational effects 
in adsorption. 

In figs. § and ^ zero-point effects due to the frustrated rotation and translation parallel 
to the surface are not taken into account since we like to concentrate here on states with 
different vibrational quantum numbers. The influence of frustrated rotation and translation 



of the H 2 molecule in contact with the surface will be discussed in a forthcoming paper j3z 

According to the principle of detailed balance, the strong enhancement of the stick- 
ing probability of vibrationally excited molecules implies a strong vibrational heating of 
molecules observed in associative desorption. Figure ^| displays the logarithm of the popu- 
lation ratio of the vibrationally first excited and the ground state in desorption versus the 
inverse surface temperature. The theoretical values were obtained by summing up over all 
final rotational states. The dashed line results from the assumption that the H 2 vibrations 
are in thermal equilibrium with the surface temperature. Indeed we find vibrational heating 
in our calculations. At T s = 700 K the ratio is 2.5 times higher than for thermal equilibrium. 
Absolute values of the vibrational heating in hydrogen desorption from Pd(100) were only 
measured for D 2 , not for H2. In the system H2/Ni(110), which has a similar sticking curve as 



H 2 /Pd (fOO) JTHj, a vibrational heating of a factor of two was found |33[, in good agreement 
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FIG. 5. Vibrational excitation in desorption. Boxes: experimental results (see text) fl8|| . 
Dashed line: vibrational population in thermal equilibrium with the surface temperature, solid 
line: 6-D calculations. 

with our results. The experimental results in Fig. |5| are only determined within a calibration 
factor [IB]. This means that only the slope of the experimental data has significance, not 
the absolute values. From this slope an apparent activation energy E a = 428 ± 30 meV 
[rj| has been deduced. Our theoretical work gives E a = 519 ± 1 meV. Considering the 
scatter in the experimental data (see Fig. ||]) and the fact that the theory does not employ 
any empirical parameters, the comparison of the experimentally and theoretically obtained 
apparent barrier is satisfactory. 

In fact, an analysis of the measured sticking probability of H2/Pd(100) ]15|| makes it 
plausible that the apparent activation energy should be close to the gas-phase frequency 
of H2. If we restrict ourselves for the sake of clarity to the two coordinates Z and d, then 
according to the principle of detailed balance the vibrational population ratio in desorption 
is given by 



Pi J e - E l kBT ° S{E,Ui = l)dE 



exp 



vib 



(4) 



Po / e~ E l k BTs s(E, Vi = 0) dE ~' r \ k B T, / 
where S(E, is the sticking probability for initial vibrational state z/j and kinetic energy E. 
Now an analysis of Eq. ^ yields that a strong lowering of the apparent activation energy from 
the gas-phase vibrational energy can only be caused by a ratio S(E, z/j = l)/S(E,b>i = 0) 
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which decreases exponentially with increasing kinetic energy. This would only be the case if 
S(E, Vi = 0) increased exponentially with increasing energy, which requires the assumption 
of a minimum barrier towards dissociative adsorption, as was done in the calculations of 



refs. However, the sticking in the system H2/Pd(100) is non-activated and initially 



decreases with increasing kinetic energy ]T5| . Such a sticking probability should not cause 
a lowering of the apparent activation energy for vibrational excitation in desorption, as 
reproduced by our calculations. Either the principle of detailed balance is not directly 
applicable for the adsorption/desorption of H 2 /Pd(100), which seems to be improbable 
from the experience of hydrogen-on- metal systems 0, or there is an inconsistency between 
the adsorption and desorption experiments. 

We further note that the apparent activation energy in our calculations depends on the 
rotational states of the molecules. If one considers, e.g., only jf = 4 rotational states, where 
jf is the final rotational quantum number, for the vibrational excitation in desorption, then 
the theoretical activation energy drops to E a = 487 ± 2 meV. This is due to the fact that 
the sticking probability for molecules in the vibrational ground state becomes small (less 
than 0.1) for ji > 4 P5| , while for vibrationally excited molecules the sticking probability is 



almost independent of the initial rotational state. Also the absolute value of the vibrational 
heating depends sensitively on the rotational quantum number; for jf = 4 the vibrational 
heating rises to a factor of 3.6. 

For D 2 /Pd(100) a vibrational over-population in desorption of v = 1 by a factor of 
nine was found ||14|| , which is much higher than our result for H 2 /Pd(100). Note that the 



interaction of D 2 with Pd(100) is given by exactly the same PES as for H 2 . At present, 
full quantum dynamical calculations for D 2 are not feasible because the number of relevant, 
energetically accessible channels is significantly higher than for H 2 due to the larger mass 
and the therefore smaller energetic level spacings. We do not expect, however, that the 
calculations for D 2 would yield a vibrational heating of a factor of nine. Again analyzing 
Eq. [|, for such a large vibrational heating the sticking probability for z/j = states of D 2 
has to be below 0.1 for all energies, much lower than for H 2 , but already the hole model, 
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that does not take into account the steering effect, yields higher values (see Fig. f|). 

The reason for this discrepancy may lie in the determination of the PES. On the other 
hand, if all experimental results were correct, the application of the principle of detailed 
balance for the adsorption/desorption of hydrogen/Pd (100) would yield a large isotope 
effect, which, for example, has not been found for hydrogen on Pd(lll) |26|| . It may also 
be that our understanding of the adsorption/desorption dynamics is still incomplete. 



III. CONCLUSIONS 

In conclusion, we reported a six- dimensional quantum dynamical study of dissociative 
adsorption on and associative desorption from H2/Pd(100). We have shown that large vi- 
brational effects in dissociative adsorption and associative desorption of H 2 /Pd(100) exist. 
They are not due to a strongly curved reaction path and a late minimum barrier to ad- 
sorption, as was previously assumed, but they are caused by the strong lowering of the H-H 
vibrational frequency during the adsorption and the multi-dimensionality of the relevant 
phase space with its broad distribution of barrier heights. Quantitative differences between 
experiment and theory and inconsistencies between adsorption and desorption experiments 
are identified which deserve further clarification. 
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